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ABSTRACT 

Pulsed thermocouple systems can be used to measure 
gas temperatures above the melting point of the 
thermocouple by various techniques of short term or 
intermittent exposure of the thermocouple operating 
at lower temperatures. This paper describes an 
approach which uses a thermocouple cooled by a 
small Jet of Inert gas. When a measurement Is to 
be made, the cooling Jet Is turned off and the 
thermocouple allowed to heat up to near Its melting 
point, at which time the cooling Is reapplied. The 
final temperature which the thermocouple should 
have attained Is then calculated by extrapolating 
an exponential curve fit to the data. Temperature 
measurements can be recorded and displayed in near 
real time by using modern high-speed computing sys- 
tems to perform these calculations. Examples of 
the technique applied to high temperature Jet en- 
gine combustor development are presented. 

INTRODUCTION 

A recent Investigation dealing with a search for 
a suitable noble-metal thermocouple to be used In 
hlgh-temperature (>1900 K) high-pressure 
(>20xl0^N/m^, >20 atm) Jet engine research combus- 
tors showed chat the Iridium leg of an Ir-40 Rh/Ir 
thermocouple Is destroyed by oxidation within sev- 
eral minutes. In conjunction with the search for 
possible thermocouple protective coatings and new 
thermocouple materials, a program was also Initi- 
ated to extend the temperature range of existing 
noble-metal thermocouples by the pulse technique. 
Various Investigators have developed systems 
to utilize pulse thermocouple techniques allowing 
an Inferential measurement of a gas temperature 
which Is higher than the melting points of the 
thermocouple materials used in the system. The 
system under study uses a Jet of inert gas to cool 
the thermocouple several hundred degrees below Its 
melting point. When a measurement Is to be made, 
the cooling Jet Is turned off and the thermocouple 
allowed to heat up to near Its melting point, at 
which time the cooling Is reapplied. The final 
temperature which the thermocouple should have 
attained Is then calculated by extrapolation using 
modern high-speed computing systems. It should 
also be noted that a secondary application of the 
present probe design with Its cooling system can 
also be utilized. This secondary mode Is for the 


case where the thermocouple materials are capable 
of operating at the final temperature, but must be 
cooled to protect them from oxidation when a meas- 
urement Is not taken. In this type of application, 
the computing portion of the system under study 
need not be used. This present report describes 
the probe design, computing techniques, and test 
and application results. 

THEORY 

Assume a thermocouple In a hot gas stream is sub- 
jected to a step Increase In temperature. The tem- 
perature of the thermocouple will Increase exponen- 
tially. The temperature of the thermocouple. In the 
absence of radiation and conduction, can be repre- 
sented by the following well known equation: 



where T is temperature, with subscript 0 as 
Initial temperature and subscript F as final tem- 
perature, T is the time constant, and t Is the 
time in seconds. Then letting Aj - Tp and 
A^ “ (Tp - Tg), equation (1) can be rewritten: 

T - Ai - Ace'"^'' (2) 

The problem Is to calculate the value of Aj when 
it cannot be directly measured (l.e., by letting 
t The calculation technique used Is to lin- 

earize equation (2) and perform a least squares best 
fit. (7) 

Equation (2) can be linearized by letting 


M -1/t (3) 

B = In(Ac) (4) 

Y 5 ln(Aj - T) (5) 

X 3 f (6) 

Then equation (2) can be rewritten In the form of 
the linear equation 

Y - MX + B (7) 
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If values for Y can be measured at various values 
of X, then the valu.a of M and B can be com- 
puted using linear regression analysis. The best 
fit Is determined by minimizing the sum of the 
squares of the deviations of the data points from 
the line. However, since the value of Aj Is 
unknown, the value of Y cannot be determined di- 
rectly. Therefore, In order to determine the best 
value of for a given set of data, the varia- 

tion of a correlation coefficient with respect to 
Ax is calculated with Y > 0 as a constraint. 

"^he beat Ax value Is the one that maximizes the 
correlation coefficient. 


the curve fit. Then the initial guess for Ax Is 
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( 11 ) 


where 2.0 was arbitrarily chosen to satisfy the con- 
dition Y > 0. Once the Initial value of Aj has 
been determined the value of the correlation coef- 
ficient can be determined from equation (8). 


These calculations have been programmed in Fortran 
on a digital computer. The flow diagram for the 
computer program Is shown In figure 1. Initially 
three values of r are calculated corresponding to: 


The correlation coefficient, which is used In mini- 
mizing the sum of the squares of the deviations of 
the data points from the line, Is given by 


Ai = Ajc 


( 12 ) 



( 8 ) 


where Y^ > 0 for all 1 and n “ number of data 
points taken along the thermocouple's exponential 
heating curve. The value of r will be in the 
range -1 < r ^ + 1 where the sign corresponds to 
the sign of the slope M. If r - 0 there is no 
correlation and If r ■ il there is a perfect fit 
with the data. 

The corresponding least squares approximation equa- 
tions for M and B are: 


n n n 
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COMPUTER PROGRAM 


(9) 


( 10 ) 


Ax - Axg + STEP (13) 

and 

Ai - Ajg + 2 (STEP) (14) 

where initially STEP ” 1000. The values of r 
are tested to determine If r has gore through a 
maximum. If not, all three values of Aj are In- 
creased by the step value (Ajg “ Axg STEP) and the 
process repeated. If a maximum has occurred be- 
tween the three values of r the steps of 1000 are 
reduced to 100 and the process repeated. The con- 
vergence is complete when the step equals one and 
the maximum value of r Is between the three calcu- 
lated values. The final temperature Is then; 


'^flnal 


Ajj, + STEP. 


(15) 


The computer program for calculating Tfxnal 
shown In table 1. The variables SMX, SMY, SMXY, 
SMXX, and SMYY are defined as follows; 


n 
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The computer program was set up to determine the 
vulue of Ax In equation (1). The procedure for 
determining Ax is to Initially guess at Aj and 
then Increase the guess by some Increment until the 
correlation coefficient r goes through a maximum. 

The data which are to be curve fit to the exponen- 
tial curve consists of a series of temperature 
readings right after the temperature step is 
applied. These readings are taken only a fraction 
of the way up the exponential. The initial guess 
for Ax is computed by determining the maximum 
value (Tjiaj() °f the temperature data point used In 


SMY - £ 
1-1 


SMXY - 

1-1 


n r, 

SMXX - X^ 
1-1 


SMVY “ ^ 
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PROBE DESCRIPTION 

A clo«e-up of the head detail of the temperature 
probe Is shown In figure 2 and an overall view of a 
bench-test set-up is shown In figure 3. A cross- 
section of the water-cooled probe body Is approxi- 
mately rectangular in shape, 1 cm thick and with a 
chord depth of 3.3 cm. The thermocouple Is a 
0.05 cm diameter Pt-13 Rh/Pt wedge shaped configu- 
ration 0.75 cm long. Coolant jet tubes are 0.16 cm 
0. D. by 0.11 cm I.D. tubes slightly projecting from 
the leading edge of the water-cooled support. Gas 
coolant flow to each tube Is controlled by Individ- 
ual solenoid valves mounted externally on the probe. 
In this manner, coolant flow Is selectively turned 
off to the Individual thermocouple which Is to be 
monitored by the data system. 

In order to minimize the delay time of the Initial 
portion of the exponential heating curve, the In- 
ternal volume of the gas coolant tube between the 
coolant jet and the solenoid valve Is minimized. 

TESTS AND RESULTS 

Tests were run In a high temperature tunnel' ' , a 
bench-test set-up (fig. 3), and two experimental 
combustors, through a gas temperature range from 
1090 to 1700 K, a Mach number range from 0.05 to 
0. 35 and at pressure levels from 1x10^ to 20x10^ 
N/m^. Coolant Jet pressure was varied to obtain 
cooling steps of 80 to 1100 K. In all of these 
tests the final gas temperature was kept below the 
melting point of the thermocouple material so that 
the thermocouple could follow the full heating 
curve. The previously described computational 
method was then applied to the Initial portion of 
the curve to derive the final temperature. The 
calculated final temperature Is then compared with 
the experimentally measured final temperature to . 
determine the curve fit error. 

The sections which follow in the TEST AND RESULTS 
include the following: The effect of variation In 

coolant jet pressure on the thermocouple Indication 
la described in the section labeled Coolant Jet 
Ef fectlvenesB . Prolonged cycling of the thermo- 
couple is discussed in the section labeled Cycling 
Tests . Also included In this section are some cal- 
culated gas temperature results obtained from ana- 
log thermocouple heating curves. (A high speed 
digitizer was not available In the combustor facil- 
ity In which the cycling tests were first per- 
formed). The Bench Tests section includes thermo- 
couple heating curve results in which the thermo- 
couple output was connected to a minicomputer. 

And, finally, tests of a prototype of a complete 
system Involving the pulsed thermocouple and compu- 
ter In an experimental combustor are described In 
the section labeled Complete System Test In 
Combustor . 

Coolant Jet Effectiveness 

Jet pressures of the coolant gas was varied during 


tunnel and combustor testing to determine the amount 
of cooling necessary to depress the thermocouple 
temperature Indications to various degrees. The 
data are presented In figure A. In these tests the 
free-stream pressure ranged from 1x10^ to 20x10^ 
N/m^. The data are plotted in terms of probe cool- 
ant pressure in excess of free-stream total pressure 
versus the cooling effectiveness expressed as a 
ratio of cooling temperature. The probe coolant 
pressure was measured just upstream of the solenoid 
valve. The results presented in figure A can be 
used to estimate the approximate coolant gas supply 
pressure necessary for a particular application. 

As an example, figure 4 Indicates that with a pres- 
sure difference of 6x10^ N/m^ the cooled thermo- 
couple would Indicate a temperature of about 1000 K 
in a 2000 K gas stream. 

A large step In cooling Is advantageous from two 
points of view. One Is that It provides a larger 
number of data points from which to extrapolate when 
calculating gas temperatures above the thermocouple 
melting point. A lower Indicated temperature, with 
the cooling gas applied, also provides more oxida- 
tion protection for the thermocouple. 

Cycling Tests 

Since a pulsed thermocouple Is Intended to be rap- 
idly heated and cooled through a large temperature 
range, for many cycles, a cycling test was performed 
In one of the experimental combustors operating at 
20x10^ N7m^. The thermocouple was cycled from 755 
to 1590 K at a rate c*" 5 seconds per cycle, for 
1 hour. After these 720 cycles, the thermocouple 
showed no apparent degradation and was, thereafter, 
subjected to further testing In the program with no 
failure. 

During the cycling tests, analog records were used 
to acquire data. One such test was made in the 
combustor at 1600 K and 20x10^ N/m^. The heating 
curves (620 to 1600 K) were repeated six times while 
combustor conditions were held steady. The resul- 
tant curves were analyzed by considering only the 
Initial 60 percent of the heating curve and extrapo- 
lating to the final temperature using the computer 
program. The curve fit error for each data set 
obtained from the heating curves was then deter- 
mined. The average of the calculated errors from 
the six data sets resulted in an error of 1.4 per- 
cent In predicting the final gas temperature. The 
maximum error for a single data et was 5.6 percent. 

Bench Tests 

The probe was set-up on a laboratory bench and oper- 
ated with heat input to the thermocouple being pro- 
vided by a propane torch (fig. 3). Water coolant 
was upplled to the probe support and laboratory 
service air supplied the pulse coolant. Jet coolant 
pressure, solenoid valve cycling time, and flame 
temperature could be varied. The bench tests were 
used as a convenient tool to acquire data for devel- 
oping the curve fitting technique. The probe was 
connected to a minicomputer to control the gas cool- 
ing and to record the thermocouple voltage by using 
digital sampling. 
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Th« btnch testa vere ut«d to sake five separate 
data sets. These data sets were used to evaluate 
the previously discussed curve fitting technique. 

A typical heating curve Is shown In figure 5 as fhe 
solid line. The data represented by the solid line 
represent one data net of 1280 sasiple points. Tine 
between each aanple point Is 2.5 ms. The computer 
program was used to curve fit Che heating curve of 
each of five data sets for 20, 40, 60, 80, and 
95 percent of the way up the teiq>erature step. The 
curve represented by 0 and □ In figure 5 is 
the exponential curve generated by a computer pro- 
gram. The O represent the portion of the curve 
that was curve fit. The D represent the ex- 
trapolated portion of the curve. 

Figures 5, 6, and 7 show typical curve fits for 95, 
60, and 20 percent of the temperature step. The 
curve fit error in determining the actual tempera- 
ture for each percentage of the exponential la 
shown in cable II. In the average error colun of 
table II, the error decreases as larger percentages 
of Che heating curve are used for curve fitting. 
This la an expected result. For adequate accuracy 
(l.e. , < 3 percent) at least 60 percent of the 
heating curve should be used for curve fitting. 

Even for percentages of 60 percent or greater, 
occasionally errors significantly greater Chan that 
for adequate accuracy will occur. Such is the case 
of the 8 percent error In Data Set 4 of table II. 
Such errors can result from tenqierature fluctua- 
tions Inherent In combustion gas streams which 
Influence the calculated exponential curve. This 
points out the shortcoming of the pulsed thermo- 
couple technique. 

Complete System Test In Combustor 

A prototype pulsed thermocouple system was Inte- 
grated into the running schedule of an experimental 
combustor program. The Jet coolant solenoid valve 
cycling operation and the thermocouple output were 
connected to and controlled by the test facility 
minicomputer. Gaseous nitrogen was used as the Jet 
coolant. The cemibustor was operated at 6.5x10^ 

N/u 2 pressure and the average cooling c'’cle temper- 
ature step was from 1130 K (coolant ».■> to 1680 K 
(coolant off). Data were also taaen for short 
periods of tine to measure the burner temperature 
fluctuations at the high end of the cycle (coolant 
off). Maxlmun peak to peak differences indicated 
by the thenaocouple were as great as 55 K (with a 
thermocouple time coratant of approximately 70 ms). 
This measurement was <rade to get an idea of the 
degree of fluctuation which would be superimposed 
on and perturb the heating curve. 

A typical heating curve obtained in the combustor 
Is shown in figure 8. The solid line again shows 
Che actual data. The O and □ show the 
calculated curve. Because of the limited speed of 
the digitizer used in the combustor facility, only 
six points on the heating curve were available for 
curve fitting. Two additional data sets were gen- 
erated from heating curves at the conditions of 
figure 8. The average of the calculated errors 
from the three data sets resulted in an error of 
2.1 percent In predicting the final gas tempera- 
ture. The maximum error for a single heating curve 
was 4.4 percent. 


The combustor tests of the prototype pulsed thermo- 
couple system were, by necessity, llsiltcd In scope 
in order to avoid interference with the main goals 
of the combustor program. The results are favorable 
even though the number of sets and data points per 
set were limited. 

CONCLUDING REMARKS 

The pulsed thermocoupl< may be applied in situations 
within the melting point of the thermocouple as well 
as when the gas stream is beyond the melting point. 
Cooling the thennorc-uple affords oxidation protec- 
tion against environments encountered in high pres- 
sure combustors. Protection is also afforded 
against "hard starts" and shut-downs in which the 
gas stream temperature may momentarily rise above 
the thermocouple melting point. 

From results presented herein, in order to use the 
pulsed thermocouple above its eieltlng point data 
must be obtained on the heating curve for about 
60 percent or more of the imposed temperature step 
in order to obtain adequate accuracy (error less 
than 3 percent). About a dozen points should be 
adequate for defining the data set from which the 
curve fit is calculated. In addition, because of 
temperature fluctuations encountered In combustion 
gas streams, the average of several heating curve 
data sets is required. The exponential curve fit- 
ting technique presented herein is designed to de- 
crease the magnitude of the effects of the gas tem- 
perature fluctuations superimposed on the thermo- 
couple heating curve. However, even using this 
technique, some of the experiments indicated that 
occasionally one set out of repeated sets of data 
would show a curve fit error much greater than the 
average of the remaining data sets. These larger 
errors can occur when Che period of temperature 
fluctuation of the combustor is of the order of the 
time constant of the pulsed thermocouple, thus, 
skewing the exponential heating curve to a high 
degree. A statistical analysis of the variation in 
curve fit errors for repeated sets of data indicate 
a very low probability for these occasional large 
errors, and suggest that they may be treated as 
"wild" points. In the present investigation such 
large curve fit errors were not discarded in the 
averaging and are presented as maximum curve fit 
errors . 

Although no radiation corrections were applied to 
the heating curves or final thermocouple indica- 
tions, such corrections are required in actual com- 
buator applications in order Co calculate combustor 
performance. For ti ■ case where the gas temperature 
is above the themfficouple melting point, the radia- 
tion correction would vary according to the indica- 
ted thermocouple output position on the heating 
curve . 

From the above considerations it is obvious that 
modern computers and associated techniques are re- 
quired in order for the pulsed thermocouple to be 
practical when applied to any fluctuating tempera- 
ture envlrvnTient such as combustor. Furthermore, 
these computers provide the potential for obtaining 
real-time temperature measurements when using pulsed 
thermocoL pies during combustor testing. 
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TABU t. - caffom PBOiauai 



C SOURCE. CUR VP IT 

c 

DOUBLE PRECISION NN,SMX,SMY,SHXY«SMXX/SHYY 
DOUBLE PRECISION THAX,AI,Ain,Y,NUM,DEN,STEP 
DOUBLE PRECISION R(3},B,M«Taim),TIME(100) 
REAL N1 

INTEOER N, l,J,FLAr, 

C N1 - NUMBER OP DATA POINTS 
C N - NUMBER OP DATA POINTS 

C NN « NUMBER Or DATA POINTS 

10 READ (6,20) N 

20 POmtAT (14) 

Nl • FLOAT (N) 

NN > DBLE (Nl) 

00 40 I « 1,N 
READ (6,50) T(l) 

40 READ (6,50) TIME(I) 

50 FORMAT (PS. 3) 

C 

C DETERMINE TMAX 

C 

TMAX - 0.0 
DO 60 I « 1,N 

IP (T(I).GT.TMAX) TMAX - T(l) 

60 COUTINUP. 

C 

Air. - TMAX*2.0 
STEP - 1000.0 
70 Al ■ Air. 

IP (AID. GT. 6000.0) GO TO 120 
C 

00 90 J • 1,3 
SHX - 0.0 
SMY - 0,0 
SMXX • 0.0 
SMYY ■ 0.0 
SMXY ■ 0.0 
C 

C CALCULATE R 

C 

00 80 I - 1,N 

Y » Al-T(t) 

Y ■ DLOG (Y) 

SMX • SMX»TIME(I) 

SMY • 3HY*Y 
SMXY - SMXY*TIHE(I )*Y 
SMXX - SMXX*TIHE(I )«*2.0 
80 SMYY - SHYY*Y»«2.0 
C NUM AND DEN ARE DUMMEY VARIABLES 
NUM • NN»SMXY«SMX*SMY 

DEtl ■ (NN«SMXX-SMX**2.0)»(NN*SMYY-SMY*«2.fl) 
DEN ■ DSQRT(DEN) 

R(J> - NUM/DEN 
R(J) •• DABS (R(J>) 

90 Al ■ AI*STEP 
C FLAG • INTEGER INDEX 
PLAG • -1 

IF <R(1).LE.R<2)) PLAG • FLAG+1 
IF (R(2).LE.R(3)) FLAG ■ eLAG*l 

IF CFLAG.EQ.l) GO TO 140 

IF (FLAG.EQ.O) GO TO 110 

IF (STEP.UE.l) CO TO 100 

Z - TMAX- 2.0 

IF (Al.LE.(TMAX-TZERO)) GO TO 120 
100 STEP - STEPyiO.O 
GO TO 70 

110 IF (STEP.NE.l) GO TO 100 
GO TO ISO 
120 WRITE (6,130) 

130 FORMAT (18H DOES NOT CONVERGE) 

STOP 

140 Air. • AIG*STEP 
GO TO 70 
C 

150 WRITE (6,170) AIG 

170 FORMAT (2lH FINAL TEMPERATURE - ,F«.1,6H DEG. 
STOP 
END 
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Figure 3, - Laboratory bench test. 
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Figure 4 - Coolant jet dfectiveness in terms of cooling 
temperature step ratio and prot>e coolant jet pressure 
tPj) in excess of stmm tc^al pressure (P^). 





Figure 5. - Convarison betvveen a typical bench test experimental 
curve to the curve fit lor 95 percent d the teRX)erature step. 
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Eigure 6. - Con^rison between a typical bench test e)q)eri mental 
curve to the curve fit for 60 percent erf the ten^ature step. 






